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Abstract 



In a rectilinear dual of a planar graph vertices are represented by simple rectilinear polygons and 
edges are represented by side-contact between the corresponding polygons. A rectilinear dual is called 
a cartogram if the area of each region is equal to a pre-specified weight of the corresponding vertex. 
The complexity of a cartogram is determined by the maximum number of corners (or sides) required for 
any polygon. In a series of papers the polygonal complexity of such representations for maximal planar 
graphs has been reduced from the initial 40 to 34, then to 12 and very recently to the currently best known 
10. Here we describe a construction with 8-sided polygons, which is optimal in terms of polygonal 
complexity as 8-sided polygons are sometimes necessary. Specifically, we show how to compute the 
combinatorial structure and how to refine the representation into an area-universal rectangular layout 
in linear time. The exact cartogram can be computed from the area-universal rectangular layout with 
numerical iteration, or can be approximated with a hill-climbing heuristic. 

We also describe an alternative construction for Hamiltonian maximal planar graphs, which allows us 
to directly compute the cartograms in linear time. Moreover, we prove that even for Hamiltonian graphs 
8-sided rectilinear polygons are necessary, by constructing a non-trivial lower bound example. The 
complexity of the cartograms can be reduced to 6 if the Hamiltonian path has the extra property that it is 
one-legged, as in outer-planar graphs. Thus, we have optimal representations (in terms of both polygonal 
complexity and running time) for Hamiltonian maximal planar and maximal outer-planar graphs. 
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1 Introduction 



There is a large body of work about representing planar graphs as contact graphs, i.e., graphs whose vertices 
are represented by geometrical objects with edges corresponding to two objects touching in some specified 
fashion. Typical classes of objects might be curves, line segments, or polygons. An early result is Koebe's 
1936 theorem [18] that all planar graphs can be represented by touching disks. 

In this paper, we consider contact representations of planar graphs, with vertices represented by simple 
interior-disjoint polygons and adjacencies represented by a non-trivial contact (shared boundary) between 
the corresponding polygons. We are specifically interested in the rectilinear weighted version where the 
vertices are represented by simple (axis-aligned) rectilinear polygons. This type of a representation is known 
as a rectilinear dual of the input planar graph. 

In the weighted version, the input is a planar graph G = (V, E) along with a weight function w : 
V(G) — > M + that assigns a weight to each vertex of G. A rectilinear dual is called a cartogram if the area 
of each region is equal to the pre-specified weight of the corresponding vertex. Such representations have 
practical applications in cartography [24], geography [28] and sociology [15], but also in VLSI Layout and 
floor-planning [22]. Other applications can be found in visualization of relational data, where using the 
adjacency of regions to represent edges in a graph can lead to a more compelling visualization than just 
drawing a line segment between two points [4]. 

For rectilinear duals (unweighted) and for cartograms (weighted) it is often desirable, for aesthetic, 
practical and cognitive reasons, to limit the polygonal complexity of the representation, measured by the 
number of sides (or by the number of corners). Similarly, it is also desirable to minimize the unused area 
in the representation, also known as holes in floor-planning and VLSI layouts. A given rectilinear dual is 
area-universal if it can realize a cartogram with any pre-specified set of weights for the vertices of the graph 
without disturbing the underlying adjacencies and without increasing the polygonal complexity. 

With these considerations in mind, we study the problem of constructing area-universal rectilinear duals 
and show how to compute cartograms with worst-case optimal polygonal complexity and without any holes. 

1.1 Related Work 

In our paper and in most of the other papers cited here, "planar graph" refers to an inner-triangulated planar 
graph with a simple outer-face; the former restriction is required if at most three rectilinear polygons are 
allowed to meet in a point and the latter restriction is customary to achieve that the union of all the polygons 
in the representation is a rectangle. 

Rectilinear duals (unweighted) were first studied in graph theoretic context, and then with renewed 
interest in the context of VLSI layouts and floor planning. It is known that 8 sides are sometimes necessary 
and always sufficient [13, 21, 34]. 

The case when the rectilinear polygons are restricted to rectangles has been of particular interest and 
there are several (independent) characterizations of the class of planar graphs that allows such rectangular 
duals [19, 20, 30]. A historical overview and a summary of the state of the art in the rectangle contact graphs 
literature can be found in Buchsbaum et al. [4]. 

In the above results on rectilinear duals and rectangular duals, the areas of the polygons are not consid- 
ered; that is, these results deal with the unweighted version of the problem. The weighted version dates back 
to 1934 when Raisz described rectangular cartograms [24]. Algorithms by van Kreveld and Speckmann [31] 
and Heilmann et al. [14] yield representations with touching rectangles but the adjacencies may be disturbed 
and there may also be a small distortions of the weights. Recently, Eppstein et al. [11] characterized the class 
of planar graphs that have area-universal rectangular duals. The construction of the actual cartogram, given 
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the area-universal rectilinear dual and the weight function, can be accomplished using a result by Wimer et 
al. [33], which in turn requires numerical iteration. 

The result of Eppstein et al. above is restricted to planar graphs that have rectangular duals. Going back 
to the more general rectilineal - duals, leads to a series of papers where the main goal has been to reduce the 
polygonal complexity while respecting all areas and adjacencies. De Berg et al. initially showed that 40 
sides suffice [8]. This was later improved to 34 sides [17]. In a recent paper [3] the polygonal complexity 
was reduced to 12 sides and even more recently to 10 sides [2]. 

Side contact representations of planar graphs have also been studied without the restriction to rectilinear 
polygons. In the unweighted case 6-sided polygons are sometimes necessary and always sufficient [12]. 
The constructive upper bound relies on convex 6-sided polygons. In the weighted version, where the area of 
each polygon is prescribed, examples are known for which polygons with 7 sides are necessary [29]. This 
lower bound is matched by constructive upper bound of 7 sides if holes are allowed [1], In the same paper 
it is shown that even allowing arbitrarily high polygonal complexity and holes of arbitrary size, there exist 
examples with prescribed areas which cannot be represented with convex polygons. If holes are not allowed 
then the best previously known polygonal complexity is 10, and it is achieved with rectilinear polygons [2]. 

1.2 Our Results 

Recall that the known lower bound on the polygonal complexity even for unweighted rectilinear duals is 
8 while the best known upper bound is 10. Here we present the first construction that matches the lower 
bound. Specifically, our construction produces 8-sided area-universal rectilinear duals in linear time, and 
is thus optimal in terms of polygonal complexity. The exact cartogram can be computed from the area- 
universal rectangular layout with numerical iteration, or can be approximated with a hill-climbing heuristic. 

For Hamiltonian maximal planar graphs we have an alternative construction which allows us to directly 
compute cartograms with 8-sided rectilinear polygons in linear time. Moreover, we prove that 8-sided 
rectilinear polygons are necessary by constructing a non-trivial lower bound example. If the Hamiltonian 
path has the extra property that it is one-legged, then we can reduce the polygonal complexity and realize 
cartograms with 6-sided polygons. This can be used to obtain 6-sided cartograms of maximal outer-planar 
graphs. Thus we have optimal (in terms of both polygonal complexity and running time) representations for 
Hamiltonian maximal planar and maximal outer-planar graphs. 

2 Preliminaries 

A planar graph, G = (V, E), is one that has a drawing without crossing in the plane along with an embed- 
ding, defined via the cyclic ordering of edges around each vertex. A plane graph is a fixed planar embedding 
of a planar graph. It splits the plane into connected regions called faces; the unbounded region is the outer- 
face and all other faces are called interior faces. A planar (plane) graph is maximal if no edge can be added 
to it without violating planarity. Thus each face of a maximal plane graph is a triangle. A Hamiltonian cycle 
in a graph G is a simple cycle containing all the vertices of G. A graph G is called Hamiltonian if it contains 
a Hamiltonian cycle. 

A set P of closed simple interior-disjoint polygons with an isomorphism V : V — > P is a polygonal 
contact representation of a graph if for any two vertices u, v G V the boundaries of V(u) and V(v) share 
a non-empty line-segment if and only if (u, v) is an edge. Such a representation is known as a rectilinear 
dual of the input graph if polygons in P are rectilinear. In the weighted version the input is the graph G, 
along with a weight function w : V(G) — > M + that assigns a weight to each vertex of G. A rectilineal - dual 
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is called a cartogram if the area of each polygon is equal to the pre-specified weight of the corresponding 
vertex. We define the complexity of a polygon as the number of sides it has. A common objective is to 
realize a given graph and a set of weights, using polygons with minimal complexity. 

2.1 Canonical Orders and Schnyder Realizers 

Next we briefly summarize the concepts of a "canonical order" of a planar graph [10] and that of a "Schnyder 
realizer" [27]. Let G = (V,E) be a maximal plane graph with outer vertices u, v, w in clockwise order. 
Then we can compute in linear time [7] a canonical order or shelling order of the vertices v\ = u, i>2 = v, 
Vs, . . ., v n = w, which is defined as one that meets the following criteria for every 4 < i < n. 

• The subgraph Gi-i C G induced by v\, vi, . . ., v%-\ is biconnected, and the boundary of its outer 
face is a cycle Cj_i containing the edge (u, v). 

• The vertex V{ is in the exterior face of and its neighbors in form an (at least 2-element) 
subinterval of the path Cj_i — (u, v). 

A Schnyder realizer of a maximal plane graph G is a partition of the interior edges of G into three sets 
S\, S2 and £3 of directed edges such that for each interior vertex v, the following conditions hold: 

• v has out-degree exactly one in each of Si, S2 and S3, 

• the counterclockwise order of the edges incident to v is: entering Si, leaving £2* entering S3, leaving 
Si, entering S2, leaving S3. 

Schnyder proved that any maximal plane graph has a Schnyder realizer and it can be computed in 0(n) 
time [27]. The first condition implies that Si, for i = 1,2,3 defines a tree rooted at exactly one exterior 
vertex and containing all the interior vertices such that the edges are directed towards the root. Denote by 
the parent of vertex v in tree T^. The following well-known lemma shows a profound connection 
between canonical orders and Schnyder realizers. 

Lemma 2.1 Let G be a maximal plane graph. Then the following hold. 

(a) A canonical order of the vertices of G defines a Schnyder realizer ofG, where the outgoing edges of a 
vertex v are to its first and last predecessor (where "first" is w.r.t. the clockwise order around v), and 
to its highest-numbered successor. 

(b) A Schnyder realizer with trees Si, S2, S3 defines a canonical order, which is a topological order of the 
acyclic graph S^ 1 U S^ 1 U S3, where Sj^ 1 is the tree Sk with the direction of all its edges reversed. 

3 Cartograms with 8-Sided Polygons 

In this section we show that 8-sided polygons are always sufficient and sometimes necessary for a cartogram 
of a maximal planar graph. Our algorithm for constructing 8-sided area-universal rectilinear duals has three 
main phases. In the first phase we create a contact representation of the graph G, where each vertex of G is 
represented by an upside-down T, i.e., a horizontal segment and a vertical segment. Figures l(a)-(b) show a 
maximal planar graph and its contact representation using T's, where the three ends of each T are marked 
with aiTows. In the second phase we make both the horizontal and vertical segments of each T into thin 
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polygons with A thickness for some A > 0. We then have a contact representation of G with T-shaped 
polygons as illustrated in Figure 1(c). In the third phase we remove all the unused area in the representation 
by assigning each (rectangular) hole to one of the polygons adjacent to it, as illustrated in Figure 1(d). We 
show that the resulting representation is an area-universal rectilinear dual of G with polygonal complexity 
8, as illustrated in Figure 1(e). 
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Figure 1: Construction of a rectilinear dual of a maximal planar graph with 8-sided polygons. 



3.1 Constructing Contact Representation with T's 

Our contact representation with T's is similar to the approach described by de Fraysseix et al. [9]. 

Let G be a planar graph. As mentioned earlier, we may assume that G is internally triangulated and has 
a simple outer-face. If need be, we can add two vertices (which we later choose as v\ and vi) and connect 
them to the outer-face to ensure that the graph is maximal. Now let v\, v-i, t>3, ■ ■ ., v n be a canonical order 
of the vertices in G with corresponding Schnyder trees S\, S2 and £3 rooted at v±, V2 and v n , Add to S\ the 
two edges (v%, v\) and (v n , v\) oriented towards v\ and add to £2 the edge (v n , V2) oriented towards V2- In 
what follows, we sometimes identify vertex Vi with its canonical label i. 

We assign to vertex i the T-shape Tj consisting horizontal and vertical segments hi and hi. Begin by 
placing T\ and T2 so that h\ is placed at y = 1, /12 is placed at y = 2, the topmost points of both 61 
and 62 have y-coordinate n + 1 and the leftmost point of the /12 touches b\. Next the algorithm iteratively 
constructs the contact representation by defining so that is placed at y = k and the topmost point of 6^ 
has y-coordinate <3?3 (A;) for 3 < k < n. After the k-th step of the algorithm we have a contact representation 
of Gk, and we maintain the invariant that the order of the vertical segments with non-empty parts in the 
half -plane y > k corresponds to the same circular order of the vertices along — (^1,^2)- 

Consider inserting Tf. for Vk- The neighbors Vk x , Vk 2 , ■ ■ ., Vk d of Vk in Gk-i form a subinterval of 
Ck-i — {v\,V2) and hence the corresponding vertical segments are also in the same order in the half-plane 
y > k — 1 of the representation of Gk-i- Since Vk = ( &3(vfe i ) for 1 < i < d (Lemma 2.1), the topmost 
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points of the corresponding vertical segments have y-coordinate k. As and Vk d are the parents of Vk in 
S\ and 52, the x-coordinates of and bk d define the x-coordinates of the two endpoints of hk- Let these 
coordinates be x\ and x r ; then hk is placed between the two points (xi, k), (x r , k) and bk is placed between 
the two points (x m , k), (x m , ^^(k)) with xi + 1 < x m < x r — 1. Finally for k = n, we place T n so that /j n 
touches b\ to the left and 62 to the right and the topmost point of b n has y-coordinate n + 1. 

We note here that this representation can be computed in linear time so that all coordinates are integers 
by pre-computing a topological order ir of 5f 1 U £2; then hk is the segment [-7r(<&i(£;)), 7r(<3?2(A;))] x k and 
bk is the segment 7r(k) x [fc, ^(fc)]. 

3.2 A-Fatteningof Ti's 

Let r' be the contact representation of G using T's obtained above. In this phase of the algorithm, we 
"fatten" T's so that each vertex is represented by a T-shaped polygon. We replace each horizontal segment 
hi by an axis-aligned rectangle Hi which has the same width as hi, and whose top (bottom) side is A/2 
above (below) hi, for some < A, as illustrated in Figure 2(a). Similarly, we replace each vertical segment 
bi by an axis-aligned rectangle Bi which has the same height as bi and whose left (right) side is A/2 to the 
left (right) of bi. We call this process X-fattening of T. Note that this process creates intersections of Hi 
with Bi, and B$ 2 u\ and intersection of Bi with H^ z ny We remove these intersections by replacing 

Hi by Hi — -B^j) — -B(j> 2 (j) and replacing Bi by Bi — Hi — H§ z ny The resulting layout is a contact 
representation V" of G where each vertex V{ of G is represented by the T-shaped polygon HiU Bi. 
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Figure 2: (a) A-fattening of T, and (b) subdividing a T-shaped polygon into four rectangles. 



3.3 Removing unused area 

In this step, we begin with the A-fat T-shaped polygonal layout, T", from above and assign each (rectan- 
gular) hole to a polygon adjacent to it. We start by placing an axis-aligned rectangle of minimum size that 
encloses F". This creates five new bounded holes. Note that all the holes in T" are rectangles, and each of 
them is bounded at the bottom by Hi for some vertex Uj. We assign each hole to this vertex. This assigns 
at most two holes to each vertex Vi. one hole L{ to the left of Bi, and one hole Ri to the right of Bi. Now 
for each vertex Vi, define Pi = Ti U Li U Ri. It is easy to see that Pj is an 8-sided rectilinear polygon 
since the left side of Li has the same x-coordinate as the left side of Hi and the right side of Ri has the 
same x-coordinate as the right side of H t . Thus we have a rectilinear dual, T, of G where each vertex Vi is 
represented by Pj. 
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We preferred the above description for the computation of T since it gives the reader some intuition for 
the construction. However, we note here that the coordinates of Pj could be computed directly, without 
going through T-shapes and A-fattening, using the values $>k{vi) for k = 1,2,3 and a topological order it 
of 5-f 1 U 52- To this end, we take a topological ordering of the acyclic graph 5f 1 U £2 and for each vertex 
Vi 7^ v n , we denote the index of Vi in this topological ordering by vr(z). We can fix the placement of the left 
and right side of Bi at x-coordinate Ti(i) — A/2 and ir(i) + A/2, respectively. Then for each vertex Vi, the 
8-sided rectilinear polygon Pj representing vi is defined as follows. 

• The horizontal base segment of Pi has y-coordinate i — A/2 and extends from the left to the right 
segment (defined below.) 

• The horizontal top segment of Pi has y-coordinate <3?3(i) (in case of P n , the top segment has y- 
coordinate n + A/2). 

• The vertical left segment of Pj has x-coordinate 7r(<3?i(i)) + A/2 and goes upward from the base 
segment (in case of P\ and P n , the left segment has x-coordinate 1 — A/2). 

• The vertical right segment of Pi has x-coordinate 7r(<3?2(^)) — A/2 and goes upward from the base 
segment (in case of P2 and P n , the right segment has x-coordinate n — 1 + A/2). 

• If Vi has no children in £2, then the left segment extends upward until the top segment. Otherwise let 
j be the vertex in it that comes clockwise after <3?i(i) in the order of neighbors around i. (One can see 
that vj is the child of Vi in £2 with the lowest canonical number.) In this case, the left segment extends 
upward until x-coordinate j — A/2, from which point the polygon continues rightward towards the 
left reflex corner at 7r(i) — A/2, and then upward until it meets the top segment. 

• If Vi has no children in «Si, then the right segment extends upward until the top segment. Otherwise 
let j be the vertex in it that comes counter-clockwise after in the order of neighbors around i. 
(One can see that Vj is the child of Vi in Si with the lowest canonical number.) In this case, the right 
segment extends upward until x-coordinate j — A/2, from which point the polygon continues leftward 
towards the right reflex corner at + A/2, and then upward until it meets the top segment. 

Then the union of these n polygons define the rectilinear dual V of G which is contained inside the 
rectangle [1 — A/2, n — 1 + A/2] x [1 — A/2, n + A/2]. Thus we can compute the representation in linear 
time, and by scaling the representation by a constat factor, we can make all coordinates to be integers of size 
0{n). 

3.4 Area-Universality 

A rectilinear dual T is area-universal if any assignment of areas to its polygons can be realized by a combi- 
natorially equivalent layout. Eppstein et al. [11] studied this concept for the case when all the polygons are 
rectangles and the outer-face boundary is also a rectangle (which they call a rectangular layout). They gave a 
characterization of area-universal rectangular layouts using the concept of "maximal line-segment". A line- 
segment of a layout is the union of inner edges of the layout forming a consecutive part of a straight-line. A 
line-segment that is not contained in any other line-segment is maximal. A maximal line-segment s is called 
one-sided if it is part of the side of at least one rectangular face, or in other words, if the perpendicular line 
segments that attach to its interior are all on one side of s. 
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Lemma 3.1 [11] A rectangular layout is area-universal if and only if each maximal segment in the layout 
is one-sided. 

No such characterization is known when some faces are not rectangles. Still we can use the characteri- 
zation from Lemma 3.1 to show that the rectilineal - dual obtained by the algorithm from the previous section 
is area-universal, with the following Lemma. 

Lemma 3.2 Let T be the rectilinear dual obtained by the above algorithm. Then V is area-universal. 

Proof: To show the area-universality of T, we divide all the polygons in T into a set of rectangles such 
that the resulting rectangular layout is area-universal. Specifically, we divide each polygon Pj into four 
rectangles Hi Bi, Li and Ri (as defined in the previous subsection) by adding three auxiliary segments: one 
horizontal and two vertical, as illustrated in Figure 2(b). Any horizontal segment s not on the bounding box 
belongs to some Hi (either top or bottom), and expanding it to its maximum it ends at B^r Vi ^ on the left 
and S$ 2 (j) on the right. So s is one-sided since it is a side of the rectangle Hi. Any vertical segment s not 
on the bounding box belongs to some Bi (either left or right), and expanding it to its maximum it ends at Hi 
on the bottom and -£f<j> 3 (j) on the top. So s is one-sided since it is a side of the rectangle Bi. 

Now given any assignment of areas w : V — > M + to the vertices V of G, we split w(vi) arbitrarily into 
four parts and assign the four values to its four associated rectangles. Since T* is area-universal, there exists 
a rectilinear dual of G that is combinatorially equivalent to T for which these areas are realized. Figure 1(f) 
illustrates the rectangular layout obtained from the rectilinear dual in Figure 1(e). □ 

So for any area-assignment, the rectilinear dual that we found can be turned into a combinatorially 
equivalent one that respects the area requirements. This proves our main result for maximal planar graphs. 
Omitting v\ and t>2 from the drawing still results in a cartogram where the union of all polygons is a rectan- 
gle, so the result also holds for all planar graphs that are inner triangulated and have a simple outer-face. 

Recall that the lower bound on the complexity of polygons in any rectilinear dual (and hence in any 
cartogram) is 8, as proven by Yeap and Sarrafzadeh [34]. The algorithm described in this section, along 
with this lower bound leads to our main theorem. 

Theorem 3.1 Eight-sided polygons are always sufficient and sometimes necessary for a cartogram of an 
inner triangulated planar graph with a simple outer-face. 

3.5 Feature Size and Supporting Line Set 

In addition to optimal polygonal complexity, we point out here that the 8-sided area-universal rectilinear 
layout constructed with our algorithm maximizes the feature size and reduces the number of supporting 
lines. Earlier constructions, e.g., [3, 8], often rely on "thin connectors" to maintain adjacencies, whereas 
our construction does not. Specifically, let G be a maximal planar graph with a prescribed weight function 
w : V(G) — > R + . Choose W and H such that W x H = A = YlveV(G) w ( v )- We are interested in 
cartograms within a rectangle of width W and the height H. Define w m i n = min v ^y^w{v). 

Recall that each vertex Vi is represented by the union of at most four rectangles Hi U Bi U Ri U L, L , 
with H, L and B - L non-empty. We can distribute the weight of Vi arbitrarily among them. In particular, we can 
assign zero areas to the rectangles Lj and Ri and split its weight into two equal parts to Hi and Bi. In this 
layout each original vertex is represented by rectangles Hi and Bi whose union is either a rectangle or some 
fattened T or L, and all the necessary contacts remain. Hence we can use this simplified layout to produce 
the cartogram. 



7 



The distribution of the weight of Vi in equal parts to Hi and Bi allows to bound the feature size. The 
height and width of each rectangle are bounded by H and W respectively. Its weight is at least w m i n /2. 
Therefore, the height and width of each rectangle is at least ^ — m ■ Thus the minimum feature size of 
the cartogram is at least 2m&4w H) - ^ ms * s worst-case optimal, as the polygon with the smallest weight 
might need to reach from left to right and top to bottom in the representation. We may choose W = H = 
\[A~. Then the minimum feature size is ^^s . Furthermore the rectangular layout based on the rectangles 
Hi and Bi alone yields a cartogram with at most 2n supporting lines, instead of the 3n supporting lines in 
the cartogram based on four rectangles per vertex. 

3.6 Computing the Cartogram 

The proof of Lemma 3.1 implies an algorithm for computing the final cartogram. Splitting the T-shaped 
polygons into four rectangles and distributing the weights on these rectangles yields an area-universal rect- 
angular dual. This combinatorial structure has to be turned into an actual cartogram, i.e., into a layout 
respecting the given weights. Wimer et al. [33] gave a formulation of the problem which combines flows 
and quadratic equations. Eppstein et al. [11] indicated that a solution can be found with a numerical iter- 
ation. Alternate methods also exist, based on non-linear programming [26], geometric programming [23], 
and convex programming [6]. Heuristic hill-climbing schemes converge much quicker and can be used in 
practice, at the expense of small errors [5, 16, 32]. 

3.7 Implementation and Experimental Results 

We implemented the entire algorithm, along with a force-directed heuristic to compute the final cartogram. 
We treat each region as a rectilinear "room" containing an amount of "air" equal to the weight assigned to 
the corresponding vertex. We then simulate the natural phenomenon of air pressure applied to the "walls", 
which correspond to the line segment borders in our layout. 

For each vertex Vi of G, the polygon Pi contains air with volume w(vi). If the area of Pi is Ai, then the 
pressure applied to each of the walls surrounding Pi is given by V{vi) = w ^y , in Section 3.4, we saw that 
the maximal segments of the layout are the two horizontal and the two vertical segment associated with each 
polygon. For each polygon, the horizontal segment other than the base is entirely inside the polygon, hence 
it feels no "pressure" on it. For each of the other three segments s for the polygon Pi, the "inward force" it 
feels is given by ^(s) = Ylv &v(s)\Pi v j)h\ ~ 'P{ v i)h- Here V(s) is the set of vertices other that Vi whose 
corresponding polygon touches the segment s and Zj (resp. L) denotes the length of s that is shared with 
Pi (resp. Pj). At each iteration, we consider the segment that feels the maximum pressure and let it move 
in the appropriate direction. The convergence of this scheme follows from [16]. Some sample input-output 
pairs are shown in Fig. 3; more examples and movies showing the gradual transformation can be found at 
www .cs.arizona. edu/ ~mjalam/ optocart. 

We ran a few simple experiments to test the heuristic for time and accuracy. In the first experiment we 
generated 5 graphs on n vertices with n in the range [10 — 50] and assigned 5 random weight distributions 
with weights in the range [10 — 100]. Next we ran the heuristic until the cartographic error dropped below 
1% and recorded the average time. All the averages were below 50 milliseconds, which confirms that good 
solutions can be found very quickly in practice; see Fig 4(a). In the second experiment we fixed the time 
allowed and tested the quality of the cartograms obtained within the time limit. Specifically, we generated 
5 graphs on n vertices with n in the range [10 — 50] and assigned 5 random weight distributions with 
weights in the range [10 — 100]. We allowed the program to run for 1 millisecond and recorded the average 
"cartographic error". Even with such a small time limit, the average cartographic error was under 2.5%; see 
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Figure 3: Input-output pairs: weights are assigned at random in the range [10,100], and the cartographic 
error in the output file is less than 1%. The colors indicate air-pressure: the greener a region is, the more it 
needs to shrink; the redder a regions is, the more it needs to grow. 

Fig. 4(b). Here, the cartographic error for a cartogram of a planar graph G = (V, E) is defined as in [31]: 
max v& v i\A(v) — w(v)\/w(v)), where w{v) denotes the weight assigned to v and A(v) denotes the area 
of the polygon representing v. All of the experiments were run on an Intel Core i3 machine with a 2.2GHz 
processor and 4GB RAM. 

4 Cartograms for Hamiltonian Graphs 

In this section we show that 8-sided polygons are always sufficient and sometimes necessary for a cartogram 
of a Hamiltonian maximal planar graph. We first give a direct linear-time construction with 8-sided regions 
without relying on numerical iteration or heuristics, as discussed in the previous section. We then prove that 
this is optimal by showing that 8 sides are necessary, with a non-trivial lower bound example. 

4.1 Sufficiency of 8-sided Polygons 

Let vi , . . . , v n be a Hamiltonian cycle of a maximal planar graph G. Consider a plane embedding of G 
with the edge (vi,v n ) on the triangular outer-face. The Hamiltonian cycle splits the plane graph G into 
two outer-planar graphs which we call the left graph Gi and right graph G r . Edges on the Hamiltonian 
cycle belong to both graphs. The naming is with respect to a planar drawing T of G in which the vertices 
v±, . . . ,v n are placed in increasing order along a vertical line, and the edges are drawn with y-monotone 
curves with leftmost edge (vi,v n ); see Figure 5(a). 

Lemma 4.1 Let G = (V, E) be a Hamiltonian maximal planar graph and let w : V — >■ M + be a weight 
function. Then a cartogram with 8-sided polygons can be computed in linear time. 
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(a) (b) 

Figure 4: Experimental results for graphs with 1 0-50 vertices. Each sample point corresponds to 25 graphs, 
(a) Plotting the average time it takes to reach cartographic error of 1%. (b) Plotting the average cartographic 
error achieved in 1 millisecond. 

Proof: Let ui, . . . , v n be a Hamiltonian cycle and F be the drawing defined above with {v\,v n ) on the 
outer-face. Suppose R is a rectangle of width W and height H where W x H = A = X^ev w(v). Each 
vertex v ; will be represented as the union of three rectangles, the left leg, the body Bi, and right leg of t> j. 
We set the width of the legs to A^ = w(vi)/(2H + W); see Figure 5(b). 

Our algorithms places vertices vi,...,v n in this order, and also reserves vertical strips for legs of all 
vertices that have earlier neighbors. More precisely, let Cj be all vertices with an edge (v j, v^) in Gi for 
which i < j < k. Similarly define IZj with respect to edges in G r . In the drawing T, Cj are those vertices 
above Vj for which the horizontal ray left from vj crosses an incident edge. 

We place vertices v±,...,Vj with the following invariant: The horizontal line through the top of Bj 
intersects, from left to right: (a) a vertical strip of width A& for each v^ G Cj, in descending order, (b) a 
non-empty part of the top of Bj, and (c) a vertical strip of width A& for each Vk G TZj, in ascending order. 

We start by placing B\ as a rectangle that spans the bottom of R. At the left and right end of the top of 
B\, we reserve vertical strips of width A& for each vertex in C\ and respectively. 

To place Bi, i > 1, first locate the vertical strips reserved for v,i in previous steps (since Vi G and 
i G TZi-i, there always are such strips, though they may have started only at the top of -Bj_i). Since vertical 
strips are in descending/ascending order, the strips for Vi are the innermost ones. Let Bi be a rectangle just 
above B^\ connecting these strips. Choose the height of Bi so large that it, together with the left and right 
leg inside the strips, has area w(vi); we will discuss soon why this height is positive. 

Finally, at the top left of the polygon of Vj we reserve a new vertical strip of width A^. for each vertex k 
that is in Ci — Ci-\. Similarly reserve strips for vertices in TZi — TZi-i- Using planarity, it is easy to see that 
vertices in Ci — must have smaller indices than vertices in and so this can be done such that the 
order required for the invariant is respected. 

Clearly this algorithm takes linear time and constructs 8-gons of the correct area. To see that it creates 
contacts for all edges, consider an edge (vi,V)~) with i < k in Gi (edges in G r are similar.) By definition 
k G Ci. If Vk G Ci — Ci-i, then we reserved a vertical strip for Vk when placing Vi. This vertical strip 
is used for the left leg of Vk, which hence touches Otherwise (v^ Ci — Ci-\) we have G 
At the time that was placed, there hence existed a vertical strip for v^. There also was a vertical strip 
for Vi G These two strips must be adjacent, because by planarity (and edge (v j, v^)) there can be no 

vertex Vj with i < j < k in So these strips create a contact between the two left legs of Vi and vu- 
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(a) 



(b) 



(d) 



(e) 



Figure 5: (a) A Hamiltonian maximal planar graph G, (b) an 8-sided polygon for vertex i, (c)-(d) illustration 
for the algorithm to construct a cartogram of G, (e) a cartogram of G with 8-sided polygons. 



We now discuss the choice of Aj = w{vi)/{2H + W). Each leg of vi has height < H and width 
Aj, hence area < HXi. Then the body Bi has area > w(vi) — 2H\ and width < W, hence height 
> (w(vi) — 2H\i) /W = Aj. It follows that Bi has positive height. Also all vertical strips fit: after placing 
vertex vi, we have a strip of width A& for each vertex Vk in Li and TZi, and these strips use width 

w(v k ) w(v k ) < 2E« fce y-K}wK) < 2(A-w(v t )) < 2A _ 2w(vj) =W _ 2X 

^ 2H+W 2H+W ~ 2H+W ~ 2H+W ~ 2H 2H+W 

Hence Bi has width > 2Aj and the polygon of v% has minimum feature size A*. □ 
This algorithm also guarantees a minimum feature size for the cartogram: mm v%& v \ = 2H+W ' wri ere 



UK, 



min„ g y w{vi). Choosing W = y2A and H = y/A/2, yields minimum feature size ^J^^/a ' 



4.2 Necessity of 8-sided Polygons 

While it was known that 8-sided rectilinear polygons are necessary for general planar graphs [25], here we 
show that 8-sided rectilinear polygons are necessary even for Hamiltonian maximal planar graphs. 

Lemma 4.2 Consider the Hamiltonian maximal planar graph G = (V,E) in Figure 6(a). Define w(j) = 
w(l) = D and w(v) = 5 for v £ V \ {j, I}, where D ^> 5. Then any cartogram of G with weight function 
w requires at least one 8-sided polygon. 

Proof: Assume for a contradiction that G admits a cartogram T with respect to w such that the polygons 
{P(v)} used in V have complexity at most 6. Observe that if {u,v,x} is some separating triangle in G, 
i.e., three mutually adjacent vertices whose removal disconnect the graph, then the region R uvw used for 
the inside of the separating triangle contains at least one reflex corner of the polygon P(u), P[v), or P(x). 
The 5-vertex set {a, c, e, g, i} in G is the union of the five separating triangles {a, c, g}, {a, c, e}, {c, e, g}, 
{a,e,i}, and {e,g,i} with disjoint interiors. Since all the polygons in F are either 4-sided or 6-sided, 
the union of the polygons for these five vertices has at most five reflex corners and hence each of the five 
separating triangles above contains the only reflex corner of the polygon for a, c, e, g, or i. In particular, the 
outer boundary of R ae i contains exactly one reflex corner from one of P{a), P{e) and P(i), and hence it is 
a rectangle, say 1234. By symmetry, we may assume that the reflex corner of P(i) is not used for R ae i- 
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Figure 6: (a) A maximal planar Hamiltonian graph with a weight function that requires at least one 8-sided 
polygon in any cartogram. (b) Illustration for the proof of Lemma 4.2. 



The 4-vertex set {a,i,k,m} is the disjoint union of three separating triangles {a,k,m}, {k, i,m}, 
{i, a, m} whose interiors are vertices l,j and n, respectively. Since the reflex corner of P(i) is not used 
for Raei, h also cannot be used for any of these separating triangles. Hence each of P(j), P(l) and P(n) 
contains exactly one reflex corner from P(a), P(k) and P(m). In particular, rectangle 1234 (which is R ae i) 
must contain the reflex corner of P(a). We also can conclude that P(j), P(l) and P(n) are all rectangles, 
since there are no additional reflex corners available to accommodate additional convex corners from P(j), 
P(l) and P{n). 

Assume the naming in Figure 6(b) is such that edge 12 belongs to P(i), edges 23 and 34 belong to P(a) 
and edge 41 belongs to P{e). By the adjacencies, P{k) must occupy corners 1 and 4 and P(n) must occupy 
corner 2, while corner 3 (which is the reflex corner of P(a)) could belong to n or I. 

Now consider the rectangles P(j) and P(l). If D is sufficiently big, then these two rectangles each 
occupy almost half of rectangle 1234. Therefore, either their a; -range or their y-range must overlap in their 
interior. Assume their y-range overlaps, the other case is similar. Which polygon should occupy the area 
that is between P(j) and P(l) horizontally? It cannot be k, because P{k) contains corners 1 and 2 and 
hence would obtain 2 reflex angles from P(j) and P(l). So it must be P(m), since n is not adjacent to j 
and I. But P{m) must also separate P(n) from both P(j) and P{1). Regardless of whether n or I occupies 
corner 3, this is not possible without two reflex vertices for m. Therefore either the areas are not respected 
or some polygon must have 8 sides. □ 

Lemma 4. 1 together with Lemma 4.2 yield the following theorem. 

Theorem 4.1 Eight-sided polygons are always sufficient and sometimes necessary for a cartogram of a 
Hamiltonian maximal planar graph. 

5 Cartograms with 6-sided Polygons 

Here we study cartograms with rectilinear 6-sided polygons. We first note that these are easily constructed 
for outer-planar graphs. Then we generalize this technique to other maximal planar Hamiltonian graphs. 
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5.1 Maximal Outer-planar Graphs 

Our algorithm from Lemma 4. 1 naturally gives drawings of maximal outer-planar graphs that use 6-sided 
polygons. Another linear-time algorithm for constructing a cartogram of a maximal outer-planar graph with 
6-sided rectilinear polygons is also described in [2], however, our construction based on Lemma 4.1 is much 
simpler. Any maximal outer-planar graph G can be made into a maximal Hamiltonian graph by duplicating 
G and gluing the copies together at the outer-face such that Gi = G = G r . (This graph has double edges, 
but the algorithm in Lemma 4. 1 can handle double edges as long as one copy is in the left and one in the 
right graph.) Create the drawing based on Lemma 4.1 with all vertices having double the weight, and cut it 
in half with a vertical line. This gives a drawing of G with 6-sided rectilinear polygons as desired. 

5.2 One-Legged Hamiltonian Cycles 

We now aim to find more maximal Hamiltonian graphs which have cartograms with 6-sided polygons. In 
a Hamiltonian cycle vi,...,v n , call vertex Vj two-legged if it has a neighbor v\ in G\ with i l < j — 1 and 
also a neighbor v\ in G r with i r < j — 1. Call a Hamiltonian cycle one-legged if none of its vertices is 
two-legged. In the construction from Lemma 4.1, the polygon of Vj obtains a reflex vertex on both sides 
only if it has a neighbor below vj-i on both sides, or in other words, if it is two-legged. Hence we have: 

Lemma 5.1 Let G = (V, E) be a maximal planar graph with a one-legged Hamiltonian cycle and let 
w : V — > R + be a weight function. Then a cartogram with Q-sided polygons can be computed in linear 
time. 

It is a natural question to characterize graphs that have such Hamiltonian cycles. Given a Hamiltonian 
cycle v\, . . . , v n we fix aplane embedding of G with outer triangle {vi,Vk, v n }. The following lemma gives 
characterization of graphs with such hamiltonian cycles. 

Lemma 5.2 Let vi, . . . ,v n be a Hamiltonian cycle in a maximal plane graph G with (vi,v n ) on the outer 
triangle. Define Wi := u n — i+i. Then the following conditions are equivalent: 

(a) The Hamiltonian cycle is one-legged. 

(b) For i = 2, . . . , n, edge (vi—i,Vi) is an outer edge of the graph Gi induced by by v\, 1)2, ■ ■ ■ , Vi (with the 
induced embedding.) 

(c) v n -\ is an outer vertex and vertex has at least two neighbors with a larger index for i = 1, . . . , n — 2. 

(d) W\, . . . , w n is a canonical ordering for G. 

(e) G admits a Schnyder realizer (S\, S2,S% ) in which w\, W2 and w n are the roots of Si, S2 and S3, 
respectively and every inner vertex is a leaf in Si or £2- 

Proof: (a) <^=^> (b): For i = 2, . . . , n we argue that (a) vertex Vi is not two-legged if and only if (b) holds 
for i. Indeed, vi) is an inner edge in Gi if and only if there are boundary edges (vi,Vj) and (vi,Vk) 

with j, k < i — 1 in Gi and G r , respectively. But then Vi is two-legged by definition. 

(b) <^=^> (c): Since v n is an outer vertex and G n = G, (b) holds for i = n if and only if v n ^i is an outer 
vertex. For i = 2, . . . ,n — lwe argue that (b) holds for i if and only if (c) holds for i — 1. Let v\, respectively 
v\, denote the third vertex in the inner facial triangle containing the edge (i^-i, Uj) in Gi, respectively G r . 
(Both triangles exist, since (uj-i, vi) is an inner edge in G.) Now Vi) is an inner edge in Gi if and 
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only if both, v\ and v\, have a smaller index than Ui_i, which in turn holds if and only if the index of every 
neighbor of different from i>j, is smaller than i — 1. 

(c) =>• (d): By (c) 102, to n } = {^n ; Un-i? ui} is the outer triangle of G, and moreover, G3, which 
is induced by v n , v n -i,v n -2, is a triangle. Hence the outer boundary of G3 is a simple cycle C3 containing 
the edge (w\, W2). In other words, the first condition of a canonical ordering is met for i = 4. Assuming (c) 
and the first condition for i = 4, . . . , n — 1, we show that the second and first condition holds for i and £ + 1, 
respectively. In the end, the second condition holds for i = n since w n is an outer vertex. 

First note that Wi is in the exterior face of Gj_i since w n lies in the exterior face and the path u>i, . . . ,w n 
is disjoint from vertices in Gi-i and the embedding is planar. By (c) Wi has at least two neighbors in Gj_i. 
If the neighbors would not form a subinterval of the path Cj_i \ (w\, W2), there would be a non-triangular 
inner face in Gi, which contains a vertex Wj with j > i in its interior. But then the path Wj, . . . ,w n , which 
is disjoint from Gi, would start and end in an interior and the exterior face of Gi, respectively. This again 
contradicts planarity. Thus the second condition of a canonical ordering is satisfied for i. Moreover Gi is 
internally triangulated, has a simple outer cycle Q containing the edge (wi,W2). In other words, the first 
condition holds for i + 1. 

(d) => (c): Since wi, . . . , w n is a canonical ordering, (wi, W2) is an outer edge. In particular, W2 = 
v n -i is an outer vertex. Clearly v\ has at least two neighbors and every neighbor has a larger index, i.e., (c) 
holds for i = 1. Moreover, by the second condition of a canonical ordering every vertex Vi = w n -i+i, 
for i = 2, ... ,n — 2, has at least two neighbors in G n _; = G \ Gi, which is the subgraph induced by 
v n ,--- ,Vi+l. 

(d) =>■ (e): Consider the Schnyder realizer (<Si, £2, £3) °f C defined by the canonical order w\ , . . . , w n 
according to Lemma 2.1. For i = 3, . . . , n — 1 the outer cycle Cj of Gj consists of the edge {w\, W2), the 
WjWi-path P\ in Si, and the u>j7x/2-pafh P 2 in 5 2 . Due to the counterclockwise order of edges in a Schnyder 
realizer, no vertex on Pi, respectively P 2 , has an incoming inner edge in Gi in S2, respectively S\. Thus 
considering only edges in Gi every outer vertex in Gi, different from w\, W2, is a leaf in S\ or S 2 . When 
in the canonical ordering vertex u^+i is attached to Gi, some vertices on Cj become inner vertices of Grf+i. 
Every inner edge in Gi + ±, which was not an edge in Gi is in S3. Thus every inner vertex in Gi is a leaf in 
either Si or S 2 . 

(e) ==> (d): Consider a canonical ordering wi, W2, ■ ■ ■ , w n of G defined by the Schnyder realizer 
(Si,S2,S3) according to Lemma 2.1. Then {wi, W2, 103} is a triangle, hence C3 consists of the edge 
(w\,W2), the tU3i(;i-path Pi in Si, and the W3^ 2 -path P2 in S2. For i = 4, . . . ,n the vertex i«j is at- 
tached to Gi-i. If Wi+i would have no edge to Wi then the outgoing edge of Wi in Si or S 2 is connected an 
inner vertex in the u>jU> 2 -path or WiWi-path, respectively. But this vertex would then have an incoming edge 
in both, Si and S 2 - a contradiction. □ 

Figure 7 shows an example of a maximal plane graph with a one-legged Hamiltonian cycle, the corre- 
sponding canonical ordering, and the Schnyder realizer. 

Once we have a one-legged Hamiltonian cycle, we can build a 6-sided cartogram via Lemma 5.1 in 
linear time. Alternately we could obtain from it a Schnyder wood, rooted such that every vertex is a leaf 
in Si or S 2 , and hence obtain a 6-sided cartogram via the algorithm in Section 3. However, we prefer the 
construction of Lemma 5. 1 due to its linear runtime. 

We know that not every Hamiltonian maximally planar graph admits a one-legged Hamiltonian cycle; 
for example, the graph in Figure 6 does not even admit a cartogram with 6-gons. However, we believe that 
some non-trivial subclasses of Hamiltonian maximally planar graphs are also one-legged Hamiltonian. In 
particular, we have the following conjecture: 
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Figure 7: A graph with a one-legged Hamiltonian cycle and the corresponding Schnyder realizer. 



Conjecture 5.1 Every 4-connected maximal planar graph has a one-legged Hamiltonian cycle. 

Note that by Lemma 5.2, the conjecture is equivalent to asking whether every 4-connected maximal 
graph has a Hamiltonian cycle such that taking the vertices in this order gives a canonical ordering. Such a 
result might be of use for other graph problems as well. 

6 Conclusion and Open Problems 

We presented a cartogram construction for maximal planar graphs with optimal polygonal complexity. For 
the precise realization of the actual cartogram this approach requires numerical iteration. Even though the 
simple heuristic works well in practice, a natural open problem is whether everything can be computed with 
an entirely combinatorial linear-time approach. 

We also presented such an entirely combinatorial linear-time construction for Hamiltonian maximal 
planar graphs and showed that the resulting 8-sided cartograms are optimal. Finally, we showed that if the 
graph admits a one-legged Hamiltonian cycle (for example outer-planar graphs), only 6 sides are needed. It 
remains to identify larger classes of planar graphs which are one-legged Hamiltonian and thus have 6-sided 
cartograms. We conjecture that 4-connected maximal planar graphs have this property. 

All of the constructions in this paper yield area-universal rectilinear duals with optimal polygonal com- 
plexity. While Eppstein et al. [11] characterized area-universal rectangular layouts, a similar characteriza- 
tion remains an open problem for general area-universal rectilinear layouts. 

For some classes of graphs the unweighted and weighted versions of the problem have the same polygo- 
nal complexity, as in the case of general planar graphs where we have shown that the tight bound of 8-sided 
for weighted graphs matches the tight bound for unweighted graphs. On the other hand, Hamiltonian max- 
imal planar graphs have a tight bound of 6 in the unweighted case, while we have shown that the tight 
bound is 8 in the weighted case. It would be interesting to study when the weighted version of the problem 
increases the polygonal complexity. 

In a similar vein, rectilinear representations are often desirable for practical and technical reasons (e.g., 
for VLSI layout or floor-planning). Sometimes, insisting on rectilinear representation increases the underly- 
ing polygonal complexity. For example, general (unweighted) planar graphs can be represented by 6-sided 
polygons (tight bound) while 8 are needed in the rectilinear case. For the weighted version, we also now 
know that 8 sided are sufficient in the rectilinear case, but can we improve this to 7 sides if we do not insist 
on rectilinear layouts? 
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